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Abstract 

A field theoretical method is developed which permits us to study the dy- 
namics of vortices in disordered environments. In particular, we obtain a 
self-consistent system of equations for disorder averaged quantities. Making 
use of a Hartree-type approximation, we calculate the current-voltage (I-V) 
characteristics in the flux flow regime. In order to probe incipient melting of 
the vortex lattice we propose an experiment where interference steps in the 
I-V characteristics are observed, which arise when a small ac field is superim- 
posed on the constant voltage. 
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Superconductors of type II allow the magnetic field to enter in quantized units of magnetic 
flux. Commonly one refers to the state with vortex lines present as the Shubnikov phase 
and in this phase with an applied transport current, an electric field may appear when the 
vortices are moving. 

In an ideal situation, the vortices are ordered in the form of a lattice. Material inhomo- 
geneities (static disorder Jl|) tend to prevent the free displacement of the vortices. On the 
other hand, there are also thermal fluctuations (dynamic disorder) which enhance the vortex 
motion. Thermal fluctuations are expected to play an important role in high-temperature 
superconductors @. 

It is appropriate to distinguish two regimes for the current carrying Shubnikov phase 
depending on whether the vortices are essentially at rest or in motion. The first case is 
of interest with regard to applications whereas in the second case, intrinsic properties of 
the vortices and their dynamics will dominate the physics. In order to support the above 
classification, we wish to draw attention to a recent experiment |j| where by the help of 
neutron scattering the two regimes can be observed. 

The theoretical analysis of these problems was mainly done considering static situa- 
tions or using perturbation theory for small fluctuations H . In the field of research on 
disordered systems, the techniques used for analyzing problems are the replica method Q 
and supersymmetry |J. These remove the so-called denominator problem in the partition 
function (generating functional) when the average over the disorder is taken. If one has to 
consider dynamical and interacting systems, the methods referred to above do not seem to 
work well. 

In this context, we recall a real-time method based on the formalism of Martin, Siggia, 
and Rose (MSR) |7||| that permits the calculation of averaged dynamical quantities. It was 
applied by several authors, e.g. in the fields of manifolds in random media M, polymer 



melts [|1(J, charge density waves PI , and vortices in type-II superconductors RI2| . 

We present here a reformulation of the MSR theory, which is based on the work of Corn- 
wall, Jackiw, and Tomboulis (CJT) |13|JT^| and which leads to a self-consistent procedure 



that comprises an expansion not only about the true values of the averaged fields, but also 
about the true values of the correlation and response functions. The main point of this 
method is a Legendre transform which yields the velocity as an independent variable in- 
stead of the force. We think that this choice is more suitable to describe the situation of an 
Abrikosov lattice moving in the flux-flow regime. On the other hand, using a rather direct 
Hartree approximation ||, it was shown that in the static limit ergodicity can be broken 
and aging effects may occur. 

We apply the CJT method to analyze the dissipation due to the motion of vortices in 
a disordered superconductor. Twenty years ago, this problem was considered by one of us 
(A.S.) |15j as well as by Larkin and Ovchinnikov |L(J in connection with interference steps 



discovered by Fiory ]17] , [18] in his experiments. With the formalism presented below, the 
influence both of disorder and of thermal noise on the vortex motion can be calculated in a 
systematic way. As the high-T c superconductors are built of layers, there is a fundamental 
interest in the behaviour of two-dimensional systems. In addition, the above mentioned 
steps are more pronounced if superconducting films are observed. Therefore, as a starting 
point, we restricted ourselves to the analysis of thin superconducting layers neglecting the 
vortex tilt. 

Let us now consider the classical equation of motion for a two-dimensional system of 
interacting particles {i} (position r l (t), mass m, harmonic interaction potential V^r 1 ^ /2) 
moving dissipatively (viscosity rj) in a random potential U(r) and driven by a force f(t) 
as well as by Langevin forces For the rest of the article we shall use the following 

conventions: time arguments will be written as index (e.g. r\), thermal averages are denoted 
by (-)t, and (-)u means the average with respect to the random potential, which is taken to 
be Gaussian. Thus, we have 

fdt'l R %4 ^rnri + nrl + ^ViM 

j 

= f t -VU(ri)+$l (1) 



with according mean values and correlations (U(r))u = and (U{r)U{r'))u = Tu(r — r') 



for disorder averages and ($1)t = as well as (£^/)t = 2r)kBT5^8 t t' for thermal averages. 

The probability for a given path of the particles can be expressed by the functional 
5-function, so that the generating functional for the response and correlation functions is 
given by 



Z([fl],[J l t ]) = JV[r\]V[R\}J 

x e iR»[l R »r+f+£-\7U(r)] e ir.J ^ ^) 

where e.g. vq is meant to represent Yli J2 a I dtrf % qf % . Since we have deliberately included 
the mass m of the particles, the Jacobian J = \6£/6r\ is a constant |Tj| and may thus be 
omitted. 

The above formulation, which is in essence the MSR theory, allows to perform the dis- 
order and thermal averages immediately. After that, one can write the following compact 
expression for the generating functional Z: 

Z([q])=J j V[x}e^' x e lS ^\ (3) 
S([x]) = (l/2)r • 3 • t - (i/2)i^ • d»d v T v {r - r 1 ) • R' u 

where we have introduced the symbols r = (R, r), q = (/, J), and the matrix 



' 2i7]k B T 1 R ^ 

v IA V 



(4) 



'•' now includes an additional summation over the matrix indices introduced above. Note 
the structural similiarity between this matrix and the inverse of a propagator in the Keldysh 
formalism [[2(J . 



In order to exploit the field theoretical formalism of Jackiw, Cornwall, and Tomboulis 
(CJT) |L3| we proceed by adding a quadratic source term (i/2)x • ^ • r to the exponent 
in (0). Taking the Legendre transform of the generating functional, that is, of the vacuum 
diagrams VF([q], [&]) = —i\nZ, we arrive at a set of equations which form the starting 
point of the CJT formalism (r means the averaged position and & represents the connected 
correlation and response functions of the particles): 



5W 



5W _ 1 
JK ~ 2 



r([t], [©]) = W([q], [J*]) - q . t - • * - ^Tr© . Jl; 
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(5) 

(6) 
(7) 



For the effective action T one can write 



r([t], [&}) = S([x\) - *-Tr \n3<5 + ^Tr^El & 

-^Tm + r 2 ([f],[0]). 



It has been proven by CJT to all orders, that T 2 can be represented as 



(8) 



r 2 ([f],[©]) = -iln J p[t] e>®~ W ^™) 



(9) 



2PI 



where S int means an expansion of S([x + x)) about x beginning from the third-order term 
in x, and '2PI' indicates that only two-particle irreducible vacuum graphs in a theory with 
propagators given by (J5 and vertices determined by S int are retained. 
Note that eq. (0) implies the Dyson equation 

.5T 2 



&- 1 = & + 3-2i 
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(10) 



If we set & = and q = (/,0), which is the physical situation, the above Dyson equation 
together with the equation of motion (H) constitues an exact self-consistent solution of the 
problem in consideration. 

For further progress we resort to a Hartree-type approximation, which means that in (|9|) 
the factor exp (iSi nt ) is replaced by 1 + iS,i nt . This allows us to perform the path integral 
exactly because only Gaussian integrals have to be evaluated. In physical situations, R = 0, 
and (5 has Keldysh structure, i.e. 



■ o g A ^ 

G R G K 
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(11) 



For a demonstration we apply the theory presented above to the simple system consisting 
of only one particle moving in two dimensions subject to a random potential ("'tilted rough 
surface" ') and to Langevin forces. Specifically, the applied force ft is taken to be constant 
in time; as a consequence, the mean particle position is r = vt. We may add that a constant 
mean velocity requires a sufficiently large driving force, so that the probability for having 
the particle trapped by a potential well is neglegible. In order to complete the presentation 
of the simple model, we assume that the random-potential correlator (in Fourier space) is 
of the form IV(Jfe) = IV(0) exp (-k 2 £ 2 ). 

Thus, in the overdamped limit the physical scales are determined by the disorder strength 
IV (0), the viscosity rj, as well as by the correlation length £. The self-consistent system of 
equation in this case reads 



+ (12) 



(rncu + irj) (to + zO) - Ef w 
1 — v ® v 
{mui + irj) (uj + z'0) — E|? a 

£ = QZGi P£ - 2i v k B m) , 

72/ 



S R = fdk ^ k elkvte _ lk{ gK_gK )k kg x k ^ 

J 47T 2 



Ef = -i / ^r^(fe)fe®fee^e-^-^ fc 



/j 



The last line represents the pinning force defined as /p := f — rjv. Note that we have / || v 
by symmetry. The indices 'L' and 'T' denote the longitudinal and transverse parts of the 
self-energy E, and we have used (a ® b)^ := a^fe^. 

We have solved these equations numerically by iteration, i.e. starting by rather arbitrary 
given self-energies and performing then repeatedly the above steps. For the velocity rjv > 
rV(0) 1,/2 /£ 2 , the iteration converges quite rapidly and yields the force- velocity relation shown 
in fig. [I]. We have found that a diagonal approximation of all matrices, e.g. Q R = G 1, 
leads to nearly the same curves. 



To check the quality of the Hartee-type approximation, we have examined the energy 
conservation. If we multiply the equation of motion (fl|) by the particle velocity, average over 
random potential and thermal forces, and subtract the unconnected part, we find 

vf P = V (([r-v] 2 ))-((rO) 

= -i V G£T + 2 v k B TGto a • (13) 

The last line explicitly yields the relation to the correlation and response functions. There- 
fore, we can test our approximation by inserting the corresponding functions calculated from 
eqs. ([12]) into the conservation relation. The agreement was found to be extremely well, and 
we believe that the relation holds analytically even after the Hartree-type approximation. 

For large velocities, fp and the self-energy £ vanish ~ 1/v. In this limit, a perturbation 
theory []T1] is permissible. On the other hand, for smaller velocities r/v < IV (O) 1 / 2 /^ 2 the 



fluctuations grow and finally become larger than the mean velocity. In this regime, the 
importance of the self-energy £ comes into play. 

After this preparatory discussion of a viscous particle moving on a rough surface, we will 
now study the dynamics of vortex lattices. The equation of motion ([]]) for this system can 
be cast in the form 

mf-it + Wit + J2 Dw rvt = ft~VU(l + r lt ) + & t , (14) 
i'^i 



where we have assumed harmonic vortex interactions |[21|| , and / denotes the equilibrium 
lattice sites. The dissipation is supposed to originate from a Bardeen-Stephen mechanism 
|[22|| . In terms of electrodynamical quantities the viscosity is related to the conductivity by 
7] = B(f) <7f where 0o = 7r^/e is the flux quantum. The velocity and force are related to 
voltage and current by E = B x v and f l = 2t x 4>o (Lorentz force). The mass of the 
vortices is small and therefore, it may be omitted at the last step of the calculations. 

We shall analyze the flux-flow regime where for constant Lorentz force / the mean vor- 
tex positions are supposed to be / + vt. One can go through the above sketched steps. 
In particular, we have done numerical calculations for small lattices and have been able to 



confirm the energy conservation. The expression used for the mean vortex position implies 
that the vortices experience only small deviations from the equilibrium lattice positions. 
This assumption should be reasonable for sufficiently large velocities and is confirmed ex- 
perimentally in a recent letter 0. 

Of special interest is a situation where in addition to the dc electric field one has a small 



ac voltage of frequency f2 [|17| , |Tql , so that we have ((r)) = v + wilcosQt. In that case, one 



can find steps in the Idc-Vdc curves at fields E nn > = (n' /n)(Q/2ir) y 2<p B/y/3. This effect is 
based on the interference of the applied oscillations with the motion of the periodic vortex 
lattice. It is completely analogous to the mechanism responsible for the occurence of Shapiro 



steps |23| in the I-V characteristics of Josephson junctions. There, a mixing of modes occurs 
due to the nonlinear (periodic) relation of phase and current. For large junctions one finds 
interference steps for n'Ml = n2eVd c - This relation also has been pointed out by Martinoli 



24| , who considered periodic pinning potentials. 

In the case of disorder, the mixing of modes originates from the nonlinearity of the 
random pinning potential. Thus, one can gain information about the physics of the vortex 
lattice. From width and height of the interference steps one can infer the shear modulus Cqq 
and the correlation Tu(r) of the pinning potential. This was done by Fiory [[T7[ using the 
perturbation theory of Schmid and Hauger ||15|| . With our method, we may derive improved 
expressions that are applicable at lower velocities. 

For 47r£c 66 ^> r/v > (0 o /47ri?) 1/ ' 2 r[/(O) 1//2 /^ 2 one has distinguishable steps and may ap- 
proximate the self-consistent scheme by the first nontrivial iteration. In addition, we neglect 
thermal fluctuations and assume the ac amplitude to be sufficiently small, so that one may 
restrict oneself to the consideration of the steps E n \ in an approximation, we take the 
matrix £ from the dc case. 

In order to have measurable effects of first order in w, one observes the in-phase resistivity 
p(E) [[Tj]]. Then, the height of the steps is given by 

hn = 32^mo^ Tui9lW±)Wmn (15) 
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with w mn = m 2 + mn + n 2 and gi = 27r(2i?/v / 30o) 1 ^ 2 - If we expand p(E) linearly near a 
step voltage where E ~ E nl + SE n /2, we may define the halfwidth of the steps by 5E n = 
h n /{2p'{E nl )). Thus, 



^ 1 In 

i 32%c 6e pfa n n Y,m^u{giwJin)w mn 
x {1 - y^^mn exp (o4J erfc(x m „)} 



, > 4 



•'•„,„ : = — 3i_^nPj__n s ^Y u {g l m')m' 2 . (16) 



This result differs from the earlier one [jn| by fluctuation effects which appear here in a 
systematic way. For instance, we obtain {(G^ — G^, t _ fl ) = (([fit - r «'t' ~~ (( r it — r i't'))} 2 )) (X - 
\/t — t' for large \t — t'\ and oc (t — t') 2 for small time differences, whereas in |15 a linear 



time dependence has been considered. In addition, the contribution of S fi is not included 



in the analysis of [15 



Note that the structure function is given by 

S(k, u) = [ dt e- iujt+ikl+ikvt e~ ik2{G oo- G ^ ; (17) 
i J 

consequently, the broadening of the interference steps, which is governed by the vortex 
fluctuations, is directly related to the structure factor. This can be seen e.g. from eqs. 
(12|). Therefore, the width of the steps can give indications for the beginning of a melting 



transition [^,0 in the vortex lattice. 



In summary, we have presented a method that permits to study the dynamics of systems 
subjected to random potentials and forces. We have combined the formalisms of MSR and 
CJT and derived a self-consistent system of equations for the averaged particle positions 
and the correlation and response functions. The advantage of this combination lies in the 
fact that it is not only an expansion about the mean positions, as it would be the case for a 
mean field theory, but also one about the full correlation and response functions. Eventually, 
we have solved the self-consistent system of equations in a Hartree type of aproximation. 



New results have been obtained for a single particle moving on a rough surface in a viscous 
medium and also, for a lattice of vortices. 

We gratefully acknowledge the helpful assistance of U. Eckern and H. Glockler. This 
work was supported by the Deutsche Forschungsgemeinschaft. 
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FIGURES 



FIG. 1. Pinning force fp(v;T) for a single particle (dimensionless units, i.e. 
»7 = U = 1,IV(0) = 1) 
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